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We present a long-wavelength approximation to the Navier-Stokes Cahn-Hilliard 
equations to describe phase separation in thin films. The equations we derive under- 
score the coupled behaviour of free-surface variations and phase separation. Since 
we are interested in the long-time behaviour of the phase-separating fluid, we restrict 
our attention to films that do not rupture. To do this, we introduce a regularising 
Van der Waals potential. We analyse the resulting fourth-order equations by con- 
structing a solution as the limit of a Galerkin approximation, and obtain existence 
and regularity results. In our analysis, we find a nonzero lower bound for the height 
of the film, which precludes the possibility of rupture. The lower bound depends 
on the parameters of the problem, and we compare this dependence with numerical 
simulations. We find that while the theoretical lower bound is crucial to the con- 
struction of a smooth, unique solution to the PDEs, it is not sufficiently sharp to 
represent accurately the parametric dependence of the observed dips in free-surface 
height. 



I. INTRODUCTION 

Below a certain critical temperature, a well-mixed binary fluid spontaneously sep- 
arates into its component parts, forming domains of pure liquid. This process can 
be characterised by the Cahm-Hilliard equation, and numerous studies describe the 
physics and mathematics of phase separation. In this paper we study phase sepa- 
ration in a thin layer, in which the the varying free surface and concentration fields 
are coupled through a pair of nonlinear evolution equations. 

Cahn and Hilliard introduced their eponymous equation in [B] to model phase 
separation in a binary alloy. Since then, the model has been used in diverse appli- 
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cations: to describe polymeric fluids [T], fluids with interfacial tension [21], and self- 
segregating populations in biology [10J. Analysis of the Cahn-Hilliard (CH) equation 
was given by Elliott and Zheng in [13], where they prove the existence, uniqueness, 
and regularity of solutions: given sufficiently smooth initial data, solutions to the 
CH equation are bounded for almost all time in the Sobolev space W 4 ' 2 . Several au- 
thors have developed generalisations of the CH equation: a variable-mobility model 
was introduced by Elliott and Garcke [13J, while nonlocal effects were considered by 
Gajweski and Zacharias [15] . These additional features do not qualitatively change 
the phase separation, and we therefore turn to one mechanism that does: the cou- 
pling of a flow field to the Cahn-Hilliard equation [7J. In this case, the Cahn-Hilliard 
concentration equation is modified by an advection term, and the flow field is either 
prescribed or evolves according to some fluid equation. Ding and co-workers [TT] 
provide a derivation of coupled Navier-Stokes Cahn-Hilliard (NSCH) equations in 
which the velocity advects the phase-separating concentration field, while concen- 
tration gradients modify the velocity through an additional stress term in the mo- 
mentum equation. Such models have formed the basis of numerical studies of binary 
fluids [3], while other studies without this feedback term highlight different regimes 
of phase separation under flow [31 [191 123] • In this paper, the NSCH equations form 
the starting point for our asymptotic analysis. 

As in other applications involving the Navier-Stokes equations, the complexity of 
the problem is reduced when the fluid is spread thinly on a substrate, and the upper 
vertical boundary forms a free surface [25]. Then, provided vertical gradients are 
small compared to lateral gradients, a long-wavelength approximation is possible, in 
which the full equations with a moving boundary at the free surface are reduced to a 
single equation for the free-surface height. In the present case, the reduction yields 
two equations: one for the free surface, and one for the Cahn-Hilliard concentration. 
The resulting thin-film Stokes Cahn-Hilliard equations have already been introduced 
by the authors in (25], although the focus there was on control of phase separation 
and numerical simulations in three dimensions. Here we confine ourselves to the 
two-dimensional case: we derive the thin-film equations from first principles, present 
analysis of the resulting equations, and highlight the impossibility of film rupture. 

Along with the simplification of the problem that thin-film theory provides, there 
are many practical reasons for studying phase separation in thin layers. Thin poly- 
mer films are used in the fabrication of semiconductor devices, for which detailed 
knowledge of film morphology is required [17] . Other industrial applications of 
polymer films include paints and coatings, which are typically mixtures of poly- 
mers. One potential application of the thin-film Cahn-Hilliard theory is in self- 
assembly [IH1 [221 EH]- Here molecules (usually residing in a thin layer) respond 
to an energy-minimisation requirement by spontaneously forming large-scale struc- 
tures. Equations of Cahn-Hilliard type have been proposed to explain the qualita- 



tive features of self-assembly [161 122], arid knowledge of variations in the film height 
could enhance these models. Indeed in [25] the authors use the present thin-film 
Cahn-Hilliard model in three dimensions to control phase separation, a useful tool 
in applications where it is necessary for the molecules in the film to form a given 
structure. 

The analysis of thin-film equations was given great impetus by Bernis and Fried- 
man in [2]. They focus on the basic thin-film equation, 



with no-flux boundary conditions on a line segment, and smooth nonnegative initial 
conditions. For n — 1 this equation describes a thin bridge between two masses of 
fluid in a Hele-Shaw cell, for n < 3 it is used in slip models as h — > [23], while for 
n — 3 it gives the evolution of the free surface of a thin film experiencing capillary 
forces [26] . Using a decaying free-energy functional, they prove the existence of non- 
negative solutions to Eq. ([T]) for n > 1, while for n > 4, and for strictly positive initial 
conditions, the solution is unique, strictly positive, and is almost always bounded 
in the W 3,2 norm. This paper has inspired other work on the subject j5j El [20], in 
which the effect of a Van der Waals term on Eq. ([I]) is investigated. These works 
provide results concerning regularity, long-time behaviour, and film rupture in the 
presence of an attractive Van der Waals force. More relevant to the present work is 
the paper by Wieland and Garcke [32] , in which a pair of partial differential equations 
describes the coupled evolution of free-surface variations and surfactant concentra- 
tion. The authors derive the relevant equations using the long-wavelength theory, 
obtain a decaying energy functional, and prove results concerning the existence and 
nonnegativity of solutions. We shall take a similar approach in this paper. 

When the binary fluid forms a thin film on a substrate, we shall show in Sec. |TT] 
that a long-wave approximation simplifies the Navier-Stokes Cahn-Hilliard equa- 
tions, which reduce to a pair of coupled evolution equations for the free surface and 
concentration. If h(x,t) is the scaled free-surface height, and c(x,t) is the binary 
fluid concentration, then the dimensionless equations take the form 





dh dJ 
dt dx 








(2a) 



where 




1 d 2 h 
C~dx~ 2 



^MJH!) 2 }> (2b) 



C 2 d / dc 
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Here C is the capillary number, r measures the strength of coupling between the 
concentration and free-surface variations (backreaction), and C n is the scaled in- 
terfacial thickness. Additionally, a is the dimensionless, spatially-varying surface 
tension, and is the body-force potential acting on the film. In this paper we take 
= — |A|/i~ 3 , the repulsive Van der Waals potential [27]. This choice stabilises the 
film and prevents rupture. Although rupture is in itself an important feature in thin- 
film equations [HI El [26], in this paper we are interested in late-time phase separation 
and it is therefore undesirable. 

We present the asymptotic analysis that converts the NSCH equations into Eq. ^ 
and prove that the model equations possess smooth solutions that are bounded for 
almost all time in the W 4 ' 2 norm. The principal tool in this analysis is the construc- 
tion of a free-energy functional for Eq. ^ that is a decaying function of time. We 
prove that h(x,t) > for all time, which is the no-rupture condition. 

The paper is organised as follows. In Sec. [IT] we discuss the Navier-Stokes 
Cahn-Hilliard equation and the scaling laws that facilitate the passage to the long- 
wavelength equations, and we derive Eq. In Sec. Ill we analyse these equations 
by constructing a decaying free-energy functional. We prove the existence of solutions 
to Eq. ([2]) and provide regularity results. We obtain a condition on the minimum 
value of the free-surface height, and show that this is never zero. Using numerical 
studies, we discuss the dependence of the minimum free-surface height on the prob- 
lem parameters in Sec. |IV[ and compare with the analytic results. Finally, in Sec. [V] 
we present our conclusions. 



II. THE MODEL EQUATIONS 

In this section we introduce the two-dimensional Navier-Stokes Cahn-Hilliard 
(NSCH) equation set. We discuss the assumptions underlying the long-wavelength 
approximation. We enumerate the scaling rules necessary to obtain the simplified 
equations. Finally, we arrive at a set of equations that describe phase separation in 
a thin film subject to arbitrary body forces. 

The full NSCH equations describe the coupled effects of phase separation and flow 
in a binary fluid. If v is the fluid velocity and c is the concentration of the mixture, 
where c = ±1 indicates total segregation, then these fields evolve as 

dv I 

— + ^ . = V ■ T - -V0, (3a) 
at p 

oc 

+ v ■ Vc = DV 2 (c 3 - c - 7 V 2 c) , (3b) 

at 

V • v = 0, (3c) 
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p r ( dvi i <9t> j ^ dc dc 



where 

1,7 ^ 9a; j 9^ J ^ dxidxj ^ ^ 

is the stress tensor, p is the fluid pressure, <fi is the body potential and p is the 
constant density. The constant v is the kinematic viscosity, v = rj/p, where r] is 
the dynamic viscosity. Additionally, /3 is a constant with units of [Energy] [Mass] _1 , 
y/j is a constant that gives the typical width of interdomain transitions, and D a 
diffusion coefficient with dimensions [Length] 2 [Time] -1 . 

If the system has a free surface in the vertical or z-direction and has infinite or 
periodic boundary conditions (BCs) in the lateral or x-direction, then the vertical 
BCs we impose are 

u = w = c z = c zzz on z = 0, (5a) 
while on the free surface z = h(x, t) they are 

hihjTij = -an, h^Ty = (5b) 
dh dh , 

W= m +U ^ (5c) 

fiidiC = 0, hidiV 2 c = 0, (5d) 

where n = {—d x h, 1)/[1 + (dxh) 2 ] 1 ^ 2 is the unit normal to the surface, t, is the unit 
vector tangent to the surface, s is the surface coordinate, a is the surface tension, 
and k is the mean curvature, 



k = V • h 



dxxh 



,21 f 



[1 + (d x h) 2 ] 

This choice of BCs guarantees the conservation of the total mass and volume, 

Mass = / c(x,t)d 2 x, Volume = / d 2 x. (6) 

JDom(i) JDom(t) 

Here Dom(t) represents the time-dependent domain of integration, owing to the 



variability of the free surface height. Note that in view of the concentration BC (5d) 



the stress BC (5b) and does not contain c(x,t) or its derivatives. 

These equations simplify considerably if the fluid forms a thin layer of mean 
thickness ho, for then the scale of lateral variations A is large compared with the 
scale of vertical variations h . Specifically, the parameter 5 = h /X is small, and after 
nondimensionalisation of Eq. ^ we expand its solution in terms of this parameter, 
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keeping only the lowest-order terms. For a review of this method and its applications, 
see [26]. For simplicity, we shall work in two dimensions, but the generalisation to 
three dimensions is easily effected [25J. 

In terms of the small parameter 8, the equations nondimensionalise as follows. 
The diffusion time scale is t = X 2 /D = h\/ {8 2 D) and we choose this to be the unit 
of time. Then the unit of horizontal velocity is uq = A/to = SD/ho so that u = 
(8D / ho) U, where variables in upper case denote dimensionless quantities. Similarly 
the vertical velocity is w = (8 2 D/ho) W. For the equations of motion to be half- 
Poiseuille at O (1) (in the absence of the backreaction) we choose p = (t)D/Kq) P and 
<p = (rjD/ho) Using these scaling rules, the dimensionless momentum equations 
are 



8Re 



dU TT dU TTr dU 

h U h w — 

dT OX dZ 



d /n ^ ,2d 2 U d 2 U 
(P + $) + 8 1 

dx 1 ' dx 2 dz 2 



2 vDdX 



2 ( 
\dX 



\dZ 



2 -i 



dc 

^DdX 



d 2 c d 2 c 1 
+ 



dX 2 dZ 2 



(7) 



5 3 Re 



dW TT dW TJ7 

+ u—— + w 



dT 



dX 
2 vD dZ 







d 


dZ 


)- 


~dZ 




' dc 




h 








dX 



(P + $) + §■ 

dc 
dZ 



l d 2 W_ 
'dX 2 ^ 

/?7 dc 
^DdZ 



,d 2 W 
~d~Z 2 

d 2 c 



d 2 c 



dX 2 + dZ 2 



where 



Re 



8_D_ 
v 



{SD/ho) h c 



0(1) 



(8) 



(9) 



The choice of ordering for the Reynolds number Re allows us to recover half-Poiseuille 
flow at O (1). We delay choosing the ordering of the dimensionless group (3"{/Dv until 
we have examined the concentration equation, which in nondimensional form is 



dc dc 
df + dX 



W 



d 2 



dx< 



(c 3 -c) 



dc 
dZ 
d 2 



dZ< 



( c 3 _ c ) _ 5 A C 2 a ^ -C 2 ^- 25 2 C 2 ° 2 



dc 



dX 4 



n dZ 4 



n dX 2 dZ 2 ' 



(10) 



where C 2 = "i/h\. By switching off the backreaction in the momentum equations 
(corresponding to fi^/Dv — > oo), we find the trivial solution to the momentum 
equations, U — W — d x (P + $) = d z (P + = 0, H — 1. The concentration 
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boundary conditions are then cz = czzz = on Z = 0, 1 which forces cz = so that 
the Cahn-Hilliard equation is simply 

or ax 2 { } 11 ax 4- 

To make the lubrication approximation consistent, we take 

6C n = C n = 6^/h = 0(1). (11) 



We now carry out a long- wavelength approximation to Eq. (10), writing U = Uq + 

0(5), W = Wq + 0(6), c = c + 6c\ + 5 2 c 2 + We examine the boundary 

conditions on c(x, t) first. They are h ■ Vc = n ■ VV 2 c = on Z = 0, H; on Z = 
these conditions are simply dzc = dzzzc = 0, while on Z = H the surface derivatives 
are determined by the relations 

n -V oc -6 2 H x dx + d z , 

h ■ VV 2 oc -5 4 H x d X xx - 5 2 H x d x d Z z + 5 2 d xx d z + d Z zz- 

Thus, the BCs on Co are simply dzCo = dzzzCo = on Z = 0, H, which forces 
c = c (X, T). Similarly, we find C\ = C\ (X,T), and 

dc 2 H x dc d 2 c 2 H x dc 

dZ =Z irdX> dZ 2 - = lfdX> foranyZG[0,tf]. 

In the same manner, we derive the results dzzzzc 2 = dzzzzC3 = 0. Using these 
facts, Eq. (10) becomes 

df + Uo dx~ 

d 2 i s \ ^2 /o 2 -,\ H x dc n/ ^ 2 d 2 H x dc ^ 2 d 4 c 4 



(cq - c ) - Clj— + (3cq - l) -ttwtp - 2C l^V2-FT^V ~ C " 



dx2 v u iv n 9X 4 v u ' h dX n 8X 2 H 8X n dZ*~ 
We now integrate this equation from Z = to H and use the boundary conditions 

<9 3 c 4 



on Z = 0, 



^ + on Z = ,ff. 



oz 3 "ox 3 A dx\ h dx dx 2 \H dx 
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After rearrangement, the concentration equation becomes 



H— + H(U ) — 



H 



2 



dX 2 
where 



-0 ~~ C ~ L'nTTTTTT — O r 



l dX< 



n H dX 



dH d 
dXdX 



l dX 2 



11 # dX 



(Uo) 



± £ U (X,Z,T)dZ 



is the vertically-averaged velocity. Introducing 

3 C 2 d 

» = C °- C °-Jidx 

the thin-film Cahn-Hilliard equation becomes 

dc dco = 1 d 

dT [ o) dX H3X 



H 



dc 
dX 



H 



dX 



(12) 



We are now able to perform the long-wavelength approximation to Eqs. ^ 
and @. At lowest order, Eq. @ is d z (P + $) = 0, since c = c (A,T), and 
hence 

P + $ = Psurf + $surf = P (A, h(x, t), T) + $ (A, t),T) . 

By introducing the backreaction strength 



0(1) 



(13) 



equation ^ becomes 



d 2 U d 



OZ 2 



9 / <9c 

dA ( surf + SUTf)+r dX [dX 



_l_ r 



<9c <9 2 c 2 
9A9Z2' 



Using dzzC2 = (H x /H) (dc /dX) this becomes 



d 2 U n 



d 



dZ 2 dX 



(P 



surf 



$surf) + 



r d 
HdX 



H 



dco 
dX 



(14) 



At lowest order, the BC (5b) reduces to 

dU _ «9£ 
~dZ ~ dX 



on Z = H, 
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which combined with Eq. ( 14 ) yields the relation 

dU dE / d r d 

~dZ-dX + (Z ~ H) \dX (Psurf + $surf) + H dX 



H 



dco 
dX 



Here S is the dimensionless, spatially- varying surface tension. Making use of the BC 
Uq = on Z = and integrating again, we obtain the result 



<9£ f d 

U Q (X, Z,T) = Z— + {\Z 2 -HZ) I — (P surf + $ surf ) + 



r d 
HdX 



H 



dco 
dX 



The vertically-averaged velocity is therefore 



(15) 



W |^ 9X 3 " y cdx2 



H l 



1 d 2 H \ r d 

$surf + 



HdX 



H 



dc 
dX 



where we used the standard Laplace- Young free-surface boundary condition to elim- 
inate the pressure, and 

Finally, by integrating the continuity equation in the Z-direction, we obtain, in a 
standard manner, an equation for free-surface variations, 



an a 

+ <tf W ) = 0. 



(17) 



Let us assemble our results, restoring the lower-case fonts and omitting ornamen- 

(18a) 



tation over the constants. The height equation (17) becomes 

dh dJ 



dt dx 



while the concentration equation (12) becomes 





where 



( c fr) _l ^ (j c ) - ® ( h^^ 
dt dx dx \ dx 



J 2^ q x ^ ydx\ C dx 2 *j hdx 



h(—\ 
\dx J 



(18b) 



;isc) 
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and 



r 2 -— 

n hdx 



dx ) 



and where we have the nondimensional constants 



C 



upD 



(18d) 



(19) 



Dv ho h a 5 2 

These are the thin-film NSCH equations. The integral quantities defined in Eq. ^ 
are manifestly conserved, while the free surface and concentration are coupled. 

We note that the relation C n = 5-^/j/ho = O (1) is the condition that the mean 
thickness of the film be much smaller than the transition layer thickness. In experi- 
ments involving the smallest film thicknesses attainable (10 -8 m) [28], this condition 
is automatically satisfied. The condition is also realised in ordinary thin films when 
external effects such as the air-fluid and fluid-substrate interactions do not prefer one 
binary fluid component or another. In this case, the vertical extent of the domains 
becomes comparable to the film thickness at late times, the thin film behaves in a 
quasi two-dimensional way, and the model equations are applicable. 

The choice of potential determines the behaviour of solutions. If interactions 
between the fluid and the substrate and air interfaces are important, the potential 
should take account of the Van der Waals forces present. A simple model potential 
is thus 

= Ah~ n , 

where A is the dimensionless Hamakar constant and typically n = 3 [26] . Here A can 
be positive or negative, with positivity indicating a net attraction between the fluid 
and the substrate and negativity indicating a net repulsion. This choice of potential 
can also have a regularising effect, preventing a singularity or rupture from occurring 



in Eq. ([18] . 

For = —\A\/h 3 (repulsive Van der Waals interactions), the system of equa- 
tions (18) has a Lyapunov functional F = F\ + F2, where 



Fi 



dx 



1 

2C 



dh 

dx 



+ 



\A\ 
2h* 



CI 



dx h 



2 



dc 
dx 



By differentiating these expressions with respect to time, we obtain the relation 



F l + F 2 



dxh 3 \ ^~ 

y ox 



1 d 2 h \A\ 
Cdx 2 + ~h? 



r d 
h dx 



h(—\ 
\dx J 



dx h 



dp 

dx 



(20) 



which is nonpositive for nonnegative h. This fact is the key to the analytic results 
of the next section. 



11 



III. EXISTENCE OF SOLUTIONS TO THE MODEL EQUATIONS 



In this section we prove that solutions to the model equations do indeed exist. We 



set C — |, r = \A\ = 1 in Eqs. (18) and focus on the resulting equation set 



dh 
di 



d_ 
dx 



f(h) 



d 3 h 
dx 3 



+ 



d_ 
dx 



1 dh 
g(h) dx 



+ d_(f(h±d_ 
dx \ g(h) dx 



(21a) 



where 



cf(h) 



d 3 h 
dx 3 



+ 



d_ 

dx 



c dh 
g{h) dx 



d \j{h) d 



+ ■0-1 g(h) 



dx \' 



d_ 
dx 



dx \ g(h) dx 
1 d 



dc 

g(h) dx \^ dx 



(21b) 



f{h) = h 3 , g(h) = h. 
The equations are defined on a periodic domain Q = [0, L\, while the initial conditions 



arc 



h(x, 0) = ho(x) > 0, c(x,0) = cq(x), h (x), c (x) e W ' (fi), 



(22) 



We shall prove that the solutions to this equation pair exist in the strong sense; 
however, we shall need the definition of weak solutions: 



A pair (h,c) is a weak solution of Eq. (21) if the following integral rela- 
tions hold: 

To r 

dt / dx ipth = 
o Jn 

T ° , f , f { -n^d 3 h 1 dh f(h) d 
dt / dxipJ -f(h) — + — — + JK 

o Jn I 



<9x 3 dx g(h) dx 



9(h) 



d f dc 



dx V dx 



and 



To 



dt J dxip t cg(h) = 
n in 

To 



<9x 3 g(/i) dx g{h) dx 



dx [dx} 



+ I dt J dxip x ^g(h) ^ 



1 d 



c — c — 



g(h) dx d.r 



dc 



, (23) 
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where To > is any time, and ip (x, t) and if) (x, t) are arbitrary differ- 
entiable test functions that are periodic on Q and vanish at t — and 
t = T . 



In a series of steps in Sees. HI A^HTG we prove this result: 



Given the initial data in Eq. (22), Eqs. (21) possess a strong solution 



endowed with the following regularity properties: 
(h,c) e L°° (0,T -W 2 ' 2 (Q))nL 2 (0,T -W 4 ' 2 (Q)) n ct-l (n x [0,T ]). 



The outline of the proof is as follows: In Sec. HI A we introduce a regularised ver 



sion of Eqs. (21), whose solution we approximate by a Galerkin sum in Sec. Ill B 



In Sec. Ill C we obtain a priori bounds on various norms of the approximate solu- 
tion. Crucially, we show that the free-surface height is always positive. This enables 
us to continue the approximate solution in the time interval [0,Tq]. In Sees. HID 



and |III E| we show that the Galerkin sum converges to a solution of the unapproxi- 
mated equations, in the appropriate limit. Finally, in Sees. Ill F and HI G we discuss 
the regularity and uniqueness properties of the solution. 



A. Regularisation of the problem 

We introduce regularised functions f E (s) and g £ (s) such that \im £ ^ f £ (s) = f(s), 
and lim e _^.o,s>o 9e{s) = g(s). For now we do not specify f £ (s), although we mention 
that a suitable choice of f £ (s) will cure the degeneracy of the fourth-order term in 
the height equation. On the other hand, we require that g e (s) have the properties: 
(i) g e (s) = s + e, for s > 0; (ii) g s (s) > 0, for s < 0; (hi) hm s _^_ 0O £? £ (s) = |e; and 
(iv) g s (s) is at least C 3 . 



From Eqs. (21), the regularised PDEs we study are 



where 



and 



h t = -J £ , x , (24a) 
(cg £ (h)) t = - (cJ e ) x - (g £ (h)v £>x ) x , (24b) 



fi £ = c 3 -c 1— (g £ (h)c x ) 

9e{h) 



Je = f £ (h)h xxx l — h x - ^j^- (g £ (h)c 2 x ) r 

g £ (h) g £ (h) 
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Equation (24b) can also be written as 



Q = ~ J £ c x - ^ { 9e {h)^ x ) x - [J £>x + g>(h)h t ] ; (25) 



this form of the concentration equation will be useful in Sec. Ill G 



B. The Galerkin approximation 

We choose a complete orthonormal basis on the interval Q, with periodic boundary 
conditions. Let us denote the basis by {0i(^)}ieN o - We consider the finite vector space 
Span{0 o , • • • ,(pn}- For convenience, let us take the <^(x)'s to be the eigenfunctions 
of the Laplacian on [0, L] with periodic boundary conditions, and corresponding 
eigenvalues —A?. Let </>o be the constant eigenf unction. We construct approximate 



solutions to the PDEs (24) as finite sums, 

n n 

h n (x,t) = y^^ n ,i(t)0i(x), c n (x,t) = y^7w,»(£)<M^)- 

i=0 i=0 

If the (smooth) initial data are given as 

oo oo 

h(x,0) = h Q (x) = y^^^ijx) > 0, c(x,0) = co(x) = ^^^(x), 

i=0 i=0 

then the initial data for the Galerkin approximation are 

n n 

h n (x,0) = h° n (x) = £t#&(*), c n (x,0) = c° n (x) = J>°<M*), 

i=0 i=0 

and the initial data of the Galerkin approximation converge strongly in the L 2 (Q) 
norm to the initial data of the unapproximated problem. Thus, there is a n 6 N 
such that h^(x) > 0, everywhere in Q, for all n > uq. Henceforth we work with 
Galerkin approximations with n > uq. 

Substitution of h n = Y^i=o Vnjfa m t° a weak form of the /i-equation yields 

^ih n ,<t> j )) = iJ e ,n,M), (26) 

where ^ f (h ) 

Je (h n , C n ) fe{hn)hn,xxx 77 7^n,:r — - \Jj E {H n ) C n>x ) x j 
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is the flux for the regularised /i-equation, and ([<p(x) , i>(x)]) is the pairing f n ipipdx. 
We recast Eq. (26) as 



dt 



{{Je,n, <t>j,x)) = ®n,j {Vn,Jn) 



where the function $„(??„, 7„) depends on 7/ n = (77^0, . . . , T]n,n) and 7„ = 
(7n,o> • • • > ln,n)i and is locally Lipschitz in its variables. This Lipschitz property arises 
from the fact that the regularised flux, evaluated at the Galerkin approximation, is 
a composition of Lipschitz continuous functions, and therefore, is itself Lipschitz 
continuous. 

Similarly, substitution of c n = Y^h=o ln,%<j>% h^o the weak form of the c- 
equation (24b) yields 



d 
dt 



{{g(hn)C n , 4>j)) = {{K e ,n , <f>j,x)) 



(27) 



where 



KsiJlniCn) c n f £ (h n )h nxxx -.h x c n . . (g £ ( h n ) c n x ) g{hn)^e,n,x i 

9e{K) 9e{K) 

Cn,Je,n 9e{hn)l^e,n,x > 

is the flux for the regularised c-equation (24b]). Rearranging gives 

{{g{h n )c n ,t , 4>j)) = ((Ke,n , 4>j,x)) ~ ((# '(h n )c n h n j , <pj)) , 

and the left-hand side can be recast in matrix form as Yli=o ^-ijin,i, where 



M, 



which is manifestly symmetric. It is positive definite because given a vector 
(£oj • • • , Cn), we have the relation 

^2d M ij£j = dxg £ i ^VnAe I (Mi) > °> for •••>£») ^ °> 

i,j Jn \ i J i,j 

which follows from the positivity of the regularised function g £ (s). We therefore have 
the following equation for J n ,j(t), 



C?7n,. 

dt 



(28) 
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Inspecting the expression for M^, g(h n ), and K EjVi , we see that r\ n and 7„ appear 
in a Lipschitz-continuous way in the expression Yll=o^ij l ((-^e,n > <&,a;))> while owing 
to the imposed smoothness of g e (s), the variables rj n , 7„ and rj n = (rj n> o, ■ ■ ■ ,Vn,n) 
appear in a Lipschitz-continuous way in the quantity ^h =q M^ ((5 , e(^n)c n ^n,i , <f>i))- 
The vector rj n can be replaced by the function $ n (r) n , j n ) and thus we obtain a 
relation 

n,j \J]m In) i 



in place of Eq. (28), where ^/ nj - (ry n ,7 n ) is Lipschitz. We therefore have a system of 
Lipschitz-continuous equations 

and thus local existence theory [TJ] guarantees a solution for the r/ n /s and 7 n) t's for 
all times t in a finite interval < t < a. This solution is, moreover, unique and 
continuous. To continue this approximate solution to the PDE problem in Eq. (21) 
up to an an arbitrary time T > 0, it is necessary to find a priori bounds on the 
approximate local-in-time solution. 



C. A priori bounds on the Galerkin approximation 

We identify the free energy 

dx [\h 2 x + G e {h)] + / dxg e (h) \\ (c 2 - l) 2 + |c; 



where G"(s) = 1/ [f e (s)g £ (s)]. Since the Galerkin approximation satisfies the weak 
form of the PDEs given in Eqs. (26) and (27), it is possible to obtain the free-energy 
decay law 



dF 
~dt 



(t) 



dxf £ (h r , 



n-Q- 



1 



KxXX + 9e{h n )h e {h n ) + g £ (h n ) 



dx 



- / dxg E (h n )fi 2 £ „ x + 



dx 



< t < a, (29) 



where f2_(t) = {x G Q\h n (x,t) < 0}. Now given the time-continuity of h n (x,t) in 
(0, cr), and the initial condition h° n (x) > (since n > no), there is a time o\ > 
such that h n (x,t) > for all x G Q and all t G (0,o"i). Therefore, = for 

t G (0,o"i), the last integral in (29) vanishes, and hence 

F[cn(x,t),h n (x,t)] < F[c n (x,0),h n (x,0)} < supF [c n (x, 0), h n (x, 0)] < oo, 
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for < t < o\. Consequently, we obtain the bound H/i^Ha — ^l, wnere < t < cri, 
and where k\ depends only on the initial conditions. We have Poincare's inequality 
for h n>x , 

1 2 / r \ 2 



\K\ 



< 



dx h„(x) 



+ 



L 

2^ 



\h II 2 

I ""n,x\\2 • 



Now J n dxh n (x,t) = Lr) n fi(t). Inspection of Eq. (26) shows that r) n)0 (t) = r] n0 (0) 
T]q. Thus, 

\\h n \\ 2 2 < L 2 \ V ° \ 2 + (J^J h^hKoo. 



Using result (Al) from Appendix |AJ we obtain the bound 

1 



I h n || oo 



\Kh + y^L\\h n:X \\ 2 = h ■ 



Additionally, the following properties hold: 

• The function h n is Holder continuous in space, with exponent | ; 

• J n dxG £ (h„) < k 4 . 

These results hold in < t < 01, and the constants ki, k 2 , k 3 , and fc 4 are independent 
of e, n, a, and a\, and in fact depend only on the functions h (x) and c (x). 

To continue the estimates to the whole interval (0, a), we need to prove that 
h n (•, > almost everywhere (a.e.). If this is true, there is a new interval [a\, a 2 ), 
<Ji < a 2 < cr, on which h n (-,t) > a.e., and we can then provide a priori bounds on 
h n (-, t) and c n (-,t) on the interval [o~i, a 2 ). It is then possible to show that h n (-., a 2 ) > 
a.e. and thus, by iteration, we extend the proof to the whole interval (0, cr), and 
find that h n (.,t) > a.e. on (0,cr). 
We have the bound 

dxG e (h n (-,t))<k 4 , (30) 



/ 

Jn 



where k^ depends only on the initial conditions, and where < t < o\. We now 
specify G e (s) in more detail. This function satisfies the condition 



G"(s) 



fe(s)g £ (s)' 



We take g s (s) to be as defined previously, and we choose a simple regularisation for 
f{s): 

Us) = g £ (s) 3 , 
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which is Lipschitz continuous. By defining 

CM ■ * 



CM 



dr GJr) 



f e (r)g e (r)' 

we obtain a function G e (s) that is positive for all s G (—00, 00), and 

, 1 



GJs) 



[s + e 



,2 ' 



s > 0. 



Using the boundedness of G e (s), and the time-continuity of h n (-,t), we apply the 
Dominated Convergence Theorem, 

lim / dx G e (h n (-,t)) = / dx lim G £ (h n (-, t)) = / dx G e (h n (•, 01)) < k^. 
Jn Jn Jn 

Similarly, since the constant k\ in the inequality ||/i n ,x||2 < ki, < t < o\ depends 

only on the initial data, we extend this last inequality to t — o~i, and thus h n (x, 01) 

is Holder continuous in space. 

In the worst-case scenario, the time o\ is the first time at which h n (x,t) touches 

down to zero, and and we therefore assume for contradiction that h n (xo,a\) = 0, 

and that h n (x, <7i) > elsewhere. Then, by Holder continuity, for any a;6fiwe have 

1 

the bound < h n (x, a\) < k\ \x — xq\ 2 , and thus 



/ dxG e {h n {-,o x y)> k -^ f 

Jn 6 J 



dx 



\x — x \ + e(2\^L + e) 



Hence, 



6 



h jn 
> -21o 



dx G e (h n (-,cri)) 
e (2\fL + e 



log 



L - x + (2VL + e]e x + \ 2\fL + e] e 



Thus, the integral J n G £ (h n (x, a\))dx can be made arbitrarily large, which contradicts 
the e- independent bound for this quantity, obtained in Eq. (30). We therefore have 



the strong condition that the set on which h n (-,ai) < is empty. Iterating the 
argument, we have the important result 



The set on which h n (-,t) < is empty, for < t < o. 
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Using the same argument, we have an estimate on the minimum value of h n (x,t), 



h min = min h n (x,t), 
xen,ie(o,cr] 



namely, 



k 2 L 

h mia + e> -kWL+\lkfL 1 " ' 



for all small positive e. Thus, 



frmin > M := -hVZ + J k 2 L + - , (31) 

a lower bound that depends only on the initial data c (x) and h (x). Note that 
this result depends on the intimate relationship between the Holder continuity of 
a function and its boundedness in the W 2 norm, a relationship that is true only in 
one dimension. Thus, generalisation of this lower bound, and by extension, long-time 
existence and uniqueness of solutions, does not necessarily hold in higher dimensions. 



Now, using Eq. (31) and the boundedness result 

/ dxg e {h n ) k (4-l) 2 + |c: 
Jn 1 



< k. 



where k$ depends only on the initial data, we obtain an a priori bound on \\c n "~ 



2 ^ ^5 

ux c 

V. 

2 



x\\2i 

2k, 



It is also possible to derive a bound on Uc^UI- We have the relation 

M' 

which gives the inequality \\c n \\f < 2||c||| + (4fc 5 /M) . Using the Holder relation 
||c||2 < |^I 3 II C I|4> we obtain a quadratic inequality in ||c|||, 



/ dx\{cl-lf< 
Jn 



., A 9 4|f2|A;5 

11 112 ~ 2 M 

with solution 

ii ii2^ mix _d__^ 

as required. From the boundedness of ||c nja ;||2 and ||c n || 2 follows the relation Halloo < 
k 6 < oo, a result that depends only on the initial conditions. Let us recapitulate 
these results: 
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• H^VxIh is uniformly bounded; 

• 1 1 h n | |oo is uniformly bounded; 

• The function h n is Holder continuous in space, with exponent |; 

• The function h n is nonzero everywhere and never decreases below a certain 
value M > 0, independent of n, e, and a. 

• || c n,a;||2 is uniformly bounded; 

• 1 1 c n | |oo is uniformly bounded; 

• The function c n is Holder continuous in space, with exponent \ . 
These results are independent of n, e and a, and hold for < t < o. 



D. Equicontinuity and convergence of the Galerkin approximation 



Using Eq. (29), we obtain the bound 



dt! J^dx^UK 



hn,xxx 



h, 



9e(h n )fe(K) 9e{.K. 



+9e(h n )d,n,x 



< F (0) , < t < a, 



a bound that is independent of n, a, and e. Since the quantity ||/ e (/i n )||oo = 
(H^nlloo + £ ) 3 is bounded above by a constant A 1 independent of n, e, and a, we 
have 



o Jn 



dt' / dxJl, 



-h 



h, 



1 



9e(h n )f £ (h n ) + g £ (h n ) 

< A ± F (0) = A 2 



and thus 



dt' ||J„J|* < Ao, 0<t<(j, 



where A 2 is independent of n, e, and a. Similarly, 

/ dt' \\ge{K)^e,nA\l < A 3, < t < a, 

Jo 



20 



and 

t 

dt' \\c n J £ , n \\l < A 2 ||c n ||^ < A 2 k%, 0<t<a, 

where A 2 and A3 are independent of n, e, and a. By rewriting the evolution equations 

as 

hn,t Je,n,xi G/evO^'){ ^e,n,xi 

where K £n = c n J £>n — g £ {h n )n £n ^ X) we see that there are uniform bounds for 
L dt' ||Je,n||2 an d J" cZt 7 H^fCe^Hl, which depend only on the initial data co(x) and 
h {x). 

Bernis and Friedman [2] proved the following claim: 

Let (fi (x,t) be a sequence of functions, each of which weakly satisfies the 
equation 

<fi,t = —Ji,x, Ji = J (<fi) ■ 
If (ft (x, •) is Holder continuous (exponent \), and if the fluxes Ji satisfy 

t 

dt' || Ji\\\ < Bi, 0<t<a, 

where B\ is a number independent of the index % and the time a, then 
there is a constant B 2 , independent of i and a, such that 

\<fi(;t 2 ) -<Pi(-,ti)\ <B 2 |t 2 -*i|*, 

for all ti and t 2 in (0, a). 

We now observe that the fluxes J £iU , and K £ ^ n satisfy the conditions of this theorem, 
and thus 

The functions h n (-,t) and c n (-,t) are Holder continuous (exponent |), for 
< t < a. 

We therefore have a uniformly bounded and equicontinuous family of functions 
no+i- We also have a recipe for constructing a uniformly bounded and 
equicontinuous approximate solution (h n (x, t), c n (x, t)), in a small interval (0,<r). 
The recipe can be iterated step-by-step, and we obtain a uniformly bounded and 
equicontinuous family of approximate solutions {(h n , c n )}^ =no+1 , on (0,T ) x Q, for 
an arbitrary time To. Then, using the Arzela-Ascoli theorem, we obtain the conver- 
gence result: 

There is a subsequence of the family {{h n , c n )}^ =no+l that converges uni- 
formly to a limit {h,c), in [0,T ] x Q. 
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We prove several facts about the pair (h, c). 



Let (h,c) be the limit of the family of functions {(h n , c n )}^l rao+1 con- 
structed in Sees. Ill A HID Then the following properties hold for this 
limit: 



1. The functions h(x,t) and c(x,t) are uniformly Holder continuous 
in space (exponent \), and uniformly Holder continuous in time 
(exponent ^); 

2. The initial condition (h, c) (x, 0) = (ho, cq) (x) holds; 

3. (h,c) satisfy the boundary conditions of the original problem (peri- 
odic boundary conditions); 

4. The derivatives {h,c) v (h,c) x , (h,c) xx , (h,c) xxx , and {h,c) xxxx are 
continuous in space and time; 

5. The function pair (h, c) satisfy the weak form of the PDEs, 




dtdx hipt + / / dtdx J £ (f x = 



f [h) (// (/')'•;). 



9e(h) g e (h) 




dtdx g £ {h)cip t + 



Qt 




dtdx K £ ip x = 0, 



K £ = cJ £ -g e (h) 



c 3 -c 



9e(h) 



(9e(h)c x ) 



where ip (x, t) and ip (x, t) are suitable test functions. 

The statements 1, 2, and 3 are obvious. Now, any pair (h n (x, t), c n (x, t)) satisfies the 
equation set 




dtdx h n d>t + 




dtdx J £ , n (j) x = 0, 



Je.n feihnjhn ,xxx , , "a , , 

g £ (h n ) g £ (h n ) 



1 t fe\hn) ( it \ 1 \ 
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dtdx g e {h n )c n ip t + / / dtdx K Ejn ip x = 0, 

' Qt q J J Qt q 



9ei.hr 



i,9c \hn) ('n.x) 



and from the boundedness of the fluxes J £j „ and K e>n in L 2 (0, T ; L 2 (Q)), it follows 
that 

(Je,n , K e>n ) - ( J £ , K £ ) , weakly in L 2 (0, T ; L 2 (fi)) , 

for a subsequence. Using the regularity theory for uniformly parabolic equations and 
the uniform Holder continuity of the (h n , c n )'s, it follows that 



The derivatives (h n , c n ) t , (h n , c n ) x , (h n , c n ) xx , (h n , c n ) xxx , and (h n , c n ) 
are uniformly convergent in any compact subset of (0,T ] x Q. 



Thus, 



Je = fe{h)h xxx - ^p^h x - (g £ {h)c 2 x ) x 



K £ = c.J E - g e (h) 



g £ (h) g £ (h) 
1 



c 3 -c 



9e(h) 



(9e(h)c x ) 



on (0, Tq] x Q, and therefore, claims 4 and 5 on p. 21 follow 



E. Convergence of regularised problem, as e — > 



The result in Sec. HID produced a solution (h e , c £ ) to the regularised problem. Due 
to the result 



h £ (x,t) > h min >M = -kiVL + \ \k\L + 



k 2 L 



2> 



>o, 



(32) 



independent of e, the argument of Sec. HID can be recycled to produce a solution 



(h, c) to the unregularised problem. This solution is constructed as the limit of a 
convergent subsequence, formally written as (h,c) = lim e ^ (h £ ,c £ ), and the results 
of the theorem in Sec. HID apply again to (h,c). The result (32) applies to h 
constructed as h = lim £ ^ h £ , and thus all the derivatives (h,c) t , (h,c) x , (h,c) xx , 
(h,c) xxx , and (h,c) xxxx are continuous on the whole space (0,T ] x Q and therefore, 
the weak solution (h, c) is in fact a strong one. 
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F. Regularity properties of the solution (h, c) 

Using a bootstrap argument, we show that the solution (h,c) belongs to the regu- 



larity classes L°° (0,T ; W 2 ' 2 (tt)) and L 2 (0,T ; W^ 2 {9)). From Sec. |mC| it follows 
immediately that 

WKh , \\c x \\2 < oo, 

with time-independent bounds. Thus, using Poincare's inequality, it follows that 
ft,, c G W 1,2 (Q) and, moreover, 

SUp \\h x \\2 , SUp 1 1 G,; || 2 < OO. 
[0,T ] [0,T ] 



From Sec. |IHD[ it follows that J and fi x belong to the regularity class 
L 2 (0, T ; L 2 (f2)), and hence J, /z x G L 2 (0, T ; The functions J, /i, and /i x 

take the form 

J h h xxx h x h h (h x c x ~\~ *2hc x c xx ^ , 
\i = c 3 - c - h' 1 (h x c x + hc xx ) , 

and 

fj, x (3c l) c x -\- h h x c x h h xx c x h h x c xx c XXXi 
respectively. We make use of the following observations: 

• The function h(x, t) is bounded from above and below, 

< h min < h(x,t) < h max < oo; 
we also have the boundedness of c(x, t), \\c\loo < oo. 

• Since fi x G L 2 (0, T ; L 2 (Q)), it follows that n G L 2 (0, T ; L 2 (Q)), by Poincare's 
inequality. 

• From this it follows that c x h x + hc xx is in the class L 2 (0, T ; L 2 (Q)). 

• Given the inequality 

1 , 



1 1 hjJjC x 1 1 2 5: ^max 1 1 1 1 oo 1 1 ^£ 1 1 2 — ^max 1 1 C x 1 1 2 



|/X|| 2 + \^\\f^ x \\2 



IVL 

we have the result h x c 2 x + hc x c xx G L 2 (0, T ; L 2 (fi)). 

Similarly, since J T ° ||/i/i/i x |||(ii: < oo, we have the bound h 2 c x + hh x c xx G 
L 2 (0,T ;L 2 (fi)). 
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Using these facts, and relegating the details to Appendix[Bj it is possible to show that 
c xx is in L 2 (0, T ; L^O)), from which follows the result h xxx , c xxx G L 1 (0, T ; L^ft)). 
These results give rise to further bounds, namely c xx G L 2 (0, T ; L 2 (ft)), and 
J T ° (it ||/i 2 Ca;||2 < oo, whence h xx G L 2 (0, T ; L 2 (ft)). Using this collection of bounds, 
we obtain 

h,ceL l (0,T 0] W 3 '\Q)) , 

and hence finally, 

h,ce L 2 (0,T ;W 3 ' 2 (Q)) . 
Thus, the solution (h, c) belongs to the following regularity class: 

(h, c) G L°° (0, T ; W 1,2 (ft)) n L 2 (0, T ; TU 3 ' 2 (ft)) n C^'S (ft x [0, T ]) . (33) 
Extra regularity is obtained by writing the equation pair as 
dh 

— + h 3 h xxxx = -3h 2 h x h xxx + ipt + ip 2 = (f (x, t) , 
dc 2 

TT7 + c xxxx = - T h x c xxx + + tp 2 = */> (a?, t) , 

where 

/ ll^ilb < sup \\c xx \\ 2 / dt|z/i|, v\ G L 2 ([0, r]) , 

JO [0,r] JO 

/ eft H^i lb < sup \\c xx \\ 2 / dt\v 2 \, V2 G L 2 ([0, r]) , 

JO [0,r] JO 

for any r G (0, To], and where ip 2 and ^2 belong to the class L 2 (0, To; L 2 (ft)). Here 
we have used the form of the concentration equation given by Eq. (25). By multi- 
plying the height and concentration equations by h xxxx and c xxxx respectively, and 
by integrating over space and time, it is readily shown that 

(h, c) G L°° (0,T ;^ 2 ' 2 (ft)), 

and hence 

(y^)GL 2 (0,T ;L 2 (ft)), 
from which follows the regularity result 

(h, c) G L°° (0, T ; IU 2 ' 2 (ft)) n L 2 (0, T ; lU 4 ' 2 (ft)) n C%'* (ft x [0, T ]) . (34) 
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G. Uniqueness of solutions 



Let us consider two solution pairs (h,c) and (h',c') and form the difference 
(Sh,6c) = (h — h',c—c!). Given the initial conditions (5c(x, 0), 5h(x, 0)) = (0,0), 
we show that (5h,5c) = (0,0) for all time, that is, that the solution we have con- 
structed is unique. We observe that that the equation for the difference 5c can be 
written in the form 

d <9 4 

dt Sc+ dx^ 5 ° = 5(P ^ X,t ^ (35) 

where 5(p(x, t) G L 2 (0, T ; L 2 (Q)), and where 5ip (5c = 0) = 0. Using semigroup 
theory [29], we find that Eq. (35) has a unique solution. Since 5c = satisfies 
Eq. (35), and since 5c (x, 0) = 0, it follows that 5c = for all times t G [0, T ]. 

It is now possible to formulate an equation for the difference 5h by subtracting 
the evolution equations of h and h! from one another, mindful that 5c = 0. We 
multiply the resulting equation by 5h xx , integrate over space, and obtain after using 
inequalities (see Appendix [Cl) 



-1 i U'2 2112 
X II oo 



2k sup ||5M2( r )< SU P ll<*M2 ( r ) / dt\\h~ l + h r 

re[0,T] ~ re[0,T] JO 

+ 4 sup H^ll^r) / dt(\\h' xxx + 4c x c xx \\l+\\h' x + h x c 2 x \\l) 
t£\o,t\ Jo 



where k and Kp are numerical constants. Using the results of Sec. Ill F it is readily 
shown that h' 1 + h' 2 c 2 x G L 2 (0,T; L°°(Q)), and that the functions h' xxx + 4c x c xx and 
^ + h x c x belong to the class L 2 (0,T; L 2 (f2)). By choosing T sufficiently small, it is 
possible to impose the inequality 



1 

2k 



I dtWh-' + ti'clWl + Kl I ^(II^+^CII^ + II^ + ^II 2 ) 

Jo Jo 
which in turn forces sup Te r 0T i ||5/ia;||i = 0; anc i hence the solution is unique. 



< 1, 



IV. PARAMETRIC DEPENDENCE OF THE HEIGHT DIP 



In this section we perform numerical simulations of the equations (18) and focus on 



one feature of the equations: the drop in the free-surface height at the boundary 



between binary fluid domains. Our results in Sec. [TTT| gave a rigorous upper bound 
for the magnitude of this height dip, as a function of the problem parameters, and 
we compare this estimate with numerical solutions. 
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We perform numerical simulations of the full equations ([18]), with initial data 
comprising a perturbation away from the unstable steady state (h,c) = (1,0). The 
free surface and concentration evolve to an equilibrium state where the salient feature 
is the formation of domains (intervals where c ~ ±1) that are separated by smooth 
transition regions, across which the free surface dips below its average value. Since 
we are interested in this characteristic asymptotic feature of phase separation, we 



shift our focus instead to the steady version from Eq. (18) obtained by setting <9 t = 
J — /4i = 0. For this reason also, we consider only the repulsive Van der Waals 
force, for which the interaction potential is given by <ft — ~~ \A\h~ 3 . We then solve 
the boundary-value problem 



1 d 2 h 

Cdx 2 



\A\C* 



1 



H^) 2 + ig) : 



(36a) 



dx 2 



c — c — 



1 dh dc 

h dx dx ' 



(36b) 



where now the domain is infinite and the boundary conditions are h (±oo) = 1, 
and (±oo) = 0. We have rescaled lengths by C n . In Fig. [I] we present numerical 
solutions exhibiting the dependence of the solutions on the parameters in Eq. (36). 
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FIG. 1: Equilibrium solutions of Eq. Q for C = C*\A\ = 1 and r = 0.1,1,10,50. In 
(a) the valley deepens with increasing r although the film never ruptures, while in (b) the 
front steepens with increasing r. (From 6 Naraigh and Thiffeault |25j.) 



As before, the free-surface height possesses peaks and valleys, where the valleys occur 
in the concentration field's transition zone. These profiles are qualitatively similar to 
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the results obtained in experiments on thin binary films [9j [30j |3T] . While the valley 
increases in depth for large r, rupture never takes place, in agreement with the theory 
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FIG. 2: The forces F cap and -FvdW for C = Cn = 1^1 = 1 an d r = 50 are shown to have 
opposite sign. 



in Sec. Ill Similar results are obtained by varying \A\, with smaller \A\ leading to 
larger dips. This behaviour is underlined by the results in Fig. [2j where we plot 
the sign and magnitude of the backreaction or capillary force F cap = —rh~ l {hc 2 x ) x 
and the Van der Waals force Fydw — |^| (h~ 3 ) x . These forces have opposite sign: 
the Van der Waals force inhibits film rupture, while the backreaction promotes film 
thinning. The thinning of binary films due to capillary or backreaction effects has 
been documented in experiments |30j . 

We are interested in the magnitude of the dip in the free-surface height, as a 
function of the problem parameters, and we plot the dependence of h min as a function 
of | A | and r in Fig. [3j We compare these results with the estimate for h m - m found in 



Sec. Ill In terms of the physical parameters of the system, this estimate is 



h m in >M = ^CLiFo + F^A]) 



e^lAI-H-Fb+FilAI) 2 
e 4C|A|-i(.Fo+^i|A|) 2 _ 1 



- 1 



where F x = \ J Q dx [h(x, 0)P ^ 0, and F = F (0) -F v The function M{\A\,C) has 
no explicit r-dependence: although F depends on r, it is possible to find initial data 
to remove this dependence. We show a representative plot of M(|j4|,C) in Fig. |lj 
Although a comparison between Fig. [3] and Fig. [4] is not exact, since the boundary 
conditions and domains are different in both cases, we see that the shape of the 
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FIG. 3: Dependence of dip magnitude h m i n on (a) the parameter \A\; (b) the parameter 
r. We have set C = C n = 1. In (a) the dip decreases in magnitude with increasing \A\, 
demonstrating the tendency of the Van der Waals force to flatten the free surface, while in 
(b) the dip increases with increasing r, which shows how the backreaction thins the film 
across transition zones. 



bound in Fig. [3] is different from that in Fig. |4} Since the bound in Fig. [3] is obtained 
from numerical simulations, and is intuitively correct, we conclude that it has the 
correct shape and that the bound of Fig. |4j while mathematically indispensable, is 
not sharp enough to be useful in determining the parametric dependence of the dip 
in free-surface height. 



V. CONCLUSIONS 



Starting from the Navier-Stokes Cahn-Hilliard equations, we have derived a pair of 
nonlinear parabolic PDEs that describe the coupled effects of phase separation and 
free-surface variations in a thin film of binary liquid. Since we are interested in the 
long-time outcome of the phase separation, we focused on liquids that experience a 
repulsive Van der Waals force, which tends to inhibit film rupture. Using physical 
intuition, we identified a decaying energy functional that facilitated analysis of the 
equations. We have shown that given sufficiently smooth initial data, solutions to the 
model equations (18) exist in a strong sense, and have certain regularity properties. 
Our proof follows the method developed in |2j] . Central to the analysis is the decaying 
free energy, and the derivation of a no-rupture condition, which prevents the film from 
touching down to zero. The no-rupture condition is valid in one spatial dimension 
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FIG. 4: A typical plot of M(\A\,C) for F = F 1 = \ and C = 1. This theoretical lower 
bound has a different shape from those in Fig. [3j which suggests that while M(|A|, C) plays 
an important role in the analysis of the model equations, it does not capture the physics 
of film thinning. 



only, and thus existence and uniqueness results will not necessarily be obtainable in 
higher dimensions. 

We carried out one-dimensional numerical simulations of the full equations and 
found that the free-surface height and concentration tend to an equilibrium state. 
The concentration forms domains; that is, extended regions where c ~ ±1. The 
domains are separated by narrow zones where the concentration smoothly transitions 
between the limiting values ±1. At the transition zones, the free surface dips below its 
mean value, a feature of binary thin-film behaviour that is observed in experiments. 
To study the magnitude of this dip as a function of the problem parameters, we 
focused on solving the equilibrium version of Eq. ( 18 ) as a boundary- value problem. 
This simplification is carried out without loss of generality, since we have shown 
that the system tends asymptotically to such a state. We have shown that the 
magnitude of the dip decreases by increasing the strength of the repulsive Van der 
Waals force, while the dip depth actually increases by increasing the strength of 
the backreaction. Thus, in the absence of the Van der Waals force, the film would 
rupture, preventing the occurrence of the phase separation so characteristic of long- 
time Cahn-Hilliard dynamics. The film-thinning tendency of the backreaction has 
been observed experimentally (9J [301 El] - Simulations involving two lateral directions 
have elsewhere been carried out by the authors [25], and the qualitative features are 
similar to those obtained here. 
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APPENDIX A 



In this appendix we give a list of nonstandard inequalities used in the paper and 
provide proofs. 

1. Let : Q C R — > R belong to the class W 1,2 {Vl). Then the following string of 
inequalities holds, 



sup |0| < ylHIi + 110,1k < ||0|| 2 + 41 110,112, 



(Al) 



Proof: Using the Fundamental Theorem of Calculus, we have 

for any distinct points x and a in Q. Since the function \(f>(x)\ is continuous on 
Q, it has a maximum value in Q, attained at the point x max . Thus, 

|0 (Smax) | = SUp |0| < ]0 (d) | + / ds — < |0 (a) | + / C?S 

fi Ja J n 



ds 

Since this is true for any a e Q, by integrating over a, we obtain the inequality 

SUp |0| < jUh + H0.ll! < ||0|| 2 + 41 110,112, 

where the last inequality follows from the monotonicity of norms. 
2. Let : Q C R -> R belong to the class lU 2 ' 1 ^). Then 

4 



110x112 < ^ll^xxlli + ylMlill^Hi- 



(A2) 



Proof: We have the identity J 4> 2 x dx — — J 00,,, true for any function with 
periodic boundary conditions. Using Holder's inequality, this becomes 



Il^ll2<ll0lloo||0 



XX 1 ■ 
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Using the relation (Al), this becomes 

1 



< 



L 



x 1 Wxx 1) 



< (tII0I|i + ^L\\<j>x\\2)\\<>,,\\i. 
which is a quadratic inequality in H^HI) with solution 



\\4>xx\\i + v/-^ll^lli + ^mun^m 

By sacrificing the sharpness of the bound, we obtain a simpler one 



Mia < £||<Mi 



1 <fe 1, 



as required. 



APPENDIX B 

In this appendix, we fill in the details missing from the discussion of the regularity 



of solutions in Sec. HI F and prove the result 

h,ceL 2 (0,T ;W 3 ' 2 (n)) . 



To begin, we notice that Sec. |III C| gives rise to the result 

\\h x \\i, \\cxh < °°' 

with time-independent bounds. Thus, using Poincare's inequality, it follows that 
h, c G W 1,2 {Q) and, moreover, 

SUp ||/ia;||2, SUp \\c x \\2 < OO. 
[0,T ] [0,3b] 



From Sec. HID it follows that J and jj x belong to the regularity class 
L 2 (0, T ; L 2 (Q)), and hence J, [i x G L 2 (0, T ; The functions J, /i, and fi x 

take the form 

J h h xxx h x h h {h x c x ~\~ t 2hc x c xx ^ , 
fi = c 3 - c - h' 1 [h x c x + hc xx ) , 

and 

fJ'x = (3c l) C x + /l /i^-Ca; /i h xx c x /i h x c xx c^^a,, (Bl) 
respectively. We make use of the following observations, 
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• The function h(x, t) is bounded from above and below, 

< h min < h(x,t) < h max < oo, 
and the boundedness of c(x,t), \\c\\oo < oo. 

• Since (i x G L 2 (0, T ; L 2 (fi)), it follows that \i G L 2 (0, T ; L 2 {VL)), by Poincare's 
inequality. 

• From this it follows that c x h x + hc xx is in the class L 2 (0, T ; L 2 (Q)). 

• Given the inequality 



1 1 hflC x 1 1 2 5: /imax 1 1 A* 1 1 oo 1 1 \ | 2 — ^max 1 1 Or 1 1 2 



we have the result /i^c 2 + hc x c xx G L 2 (0, T ; L 2 (f2)). 



Similarly, since \\hnh x \\ldt < oo, we have the bound /i 2 c x + hh x c xx G 
L 2 (0,T ;L 2 (fi)). 



Now inspection of fi shows that c xx is in L 2 (0, T ; L 1 (fi)), from which follows 
the result h xxx ,c xxx G L 1 (0, T ; By repeating the same argument, we find 

that c xx G L 2 (0, To; L 2 (f2)). We also have the result that ||/i 2 c x ||2 is almost always 
bounded. To show that h xx G L 2 (0, T ; L 2 (f2)), we take the evolution equation for 
h(x,t), multiply it by h, and integrate, obtaining 

~\~T [ dx h 2 — 3 [ dxh 2 h 2 x h xx — f dxh x J — f dxh 3 h xx , 
dt Jq Jq Jq Jn 

where Jo = h x h~ l + 2h 3 c x c xx + h 2 h x c 2 . The time-integral of the first term on the 
left-hand side of this equation is obviously bounded in time. Let us examine the 
time-integral of the second term, 



/•Jo r r-io 

/ dt dxh 2 h 2 x h xx < hlnx sup \\h x \\\ / dt 
Jo Jn [o,t ] Jo 

/•To 

< Cx SU P IIM2 / dt + ||/l X xx||l) < OO- 

[0,T ] Jo 

The third term on the left-hand side is dispatched in a similar way, so that 
Jo°dt \\h xx \\ 2 < oo. We have now shown that h,c G L 2 (0, T ; W 2 > 2 {VL)). Using 
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this result, together with the previous facts gathered together in this appendix, it is 
readily shown that h,c G L 1 (0, T ; W /3 ' 1 (fi)), and it follows that 

/i,cGL 2 (0,T ;W 3 ' 2 (n)) . 

To see this more clearly, we show that c xxx is in the above class. For example, 
consider a typical term in fj, x , 

f dt f dxh 2 x c 2 xx < sup \\h x \\l f dt {iT^c^x + \\c xxx \\i) 2 < oo, 
Jo Jn [o,t ] Jo 

which implies h x c xx G L 2 (0, To; L 2 (Q)). This bound, together with hfih x G 
L 2 (0,T ; L 2 (Q)), implies h 2 x c x G L 2 (0, T ; L 2 (fi)). Gathering all these results, we 
have 

fix j h x c xx , /i^Cj; , h xx c x G L 2 (0, T ; L 2 (f2)) 



From (Bl), it follows that c xxx is in this class as well. A similar argument holds for 
h xxx . Thus, the solution (h, c) belongs to the regularity class 

(M G L°° (0,To;W 1,2 (fi)) DL 2 (0,T ;iy 3 ' 2 (fi)) nd'i (fi x [0,T ]). 



APPENDIX C 

In this appendix, we describe in full the proof of the the uniqueness of solutions 



sketched in Sec. Ill G We consider two solution pairs (h, c) and (h', d) and form the 
difference (Sh, 8c) = {h — h',c — d). Given the initial conditions (Sc(x, 0), Sh(x, 0)) = 
(0,0), we show that (Sh,8d) = (0,0) for all time, that is, that the solution we have 
constructed is unique. We observe that that the equation for the difference 5c can 
be written in the form 

^+^c = Mx,t), (ci) 

where 5<p (x,t) G L 2 (0,T ; L 2 (f2)), and where 5(p (5c = 0) = 0. As discussed previ- 
ously, this equation has a unique solution, given smooth initial data. Since 5c = 



satisfies Eq. (CI ), and since 5c (x, 0) = 0, it follows that 5c = for all times t G [0, T ). 

It is now possible to formulate an equation for the difference 5h by subtracting the 
evolution equations of h and h' from one another, mindful that 5c = 0. We multiply 



34 



the resulting equation by Sh xx and integrate over space to obtain 



1 d 
~d~t 



2 



dx5hl+ [ dxh 3 8h 2 xxx = [ dx5hl xx 5h x (h-' + h> 2 c x ) 
Jn Jn 

- dx (h 3 - ti 3 ) (h' xxx + Ac x c xx ) 5h xxx 
Jn 

+ [ dx5h xxx [h' x (h' 1 - h'- 1 ) + h x c 2 x (h 2 - ti 2 )} . 
Jn 

Using the lower bound on h(x,t) > h min > and Young's first inequality, this 
equation is transformed into an inequality, 

\j t f n dx5h ^ + ^min J dx Sh 2 xxx < k x J^dx5h 2 xxx + j^dxSh 2 (V 1 + h' 2 c 2 x ) 2 

+ k 2 I dx5h 2 xxx + - i - [ dx (h 3 -ti 3 )\ti xxx + Ac x c xx f 
Jn 4K 2 Jn 

+ kJ dxdh 2 xxx -\- — — f dx [ti x (h- 1 - ti' 1 ) + h x c 2 x (h 2 - ti 2 )]\ 
Jn 4K 3 Jn 

where Ki, k 2 , and k 3 are arbitrary positive constants. By choosing ki+k 2 +k 3 = h 3 ^^, 
the inequality simplifies to 

2k j t J dx5h 2 x < J dx5h 2 x (h' l + ti 2 c 2 x ) 2 

+ [ dxbti 2 \(ti xxx + Ac x c xx f + (ti x + h x c 2 x ) 2 , 
Jn L J 

where k is another positive constant. We integrate over the time interval [0, T], and 
use the fact that ||5/i x ||2 (0) = to obtain 



J2 2 || 2 

X II oo 



2k sup H^ll^r) < sup H^ll^r) / dt\\h' l + ti 

re[0,T] re[0,T] " Jo 

+ sup / dt dx (ti xxx + Ac x c xx f + (ti x + h x c 2 x ) 

te[o,t] Jo Jn 1 

The Poincare inequality can be combined with the one-dimensional differential in- 
equalities discussed in Appendix A to yield the relation \\f\\oo < «p||/e||2, where / is 
some mean- zero function and Kp is an /-independent constant. We therefore arrive 
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at 



2k sup \\6h x \\ 2 2 (T) < sup \\6h x \\ 2 2 (T) [ dt + /i /2 c^||^ 

re[0,T] ~ re[0,T] Jo 

+ k 2 p sup ||5/i x ||^(r) / dt (H/4,.3 + 4:c x c xx \\ 2 2 + \\h' x + /i s c£||l) • 

re[0,Tl JO 



Using the results of Sec. |IIIF[ it is readily shown that h~ l + h' 2 c 2 x G L 2 (0, T; 

and that the functions h' xxx +Ac x c xx and h' x + h x c 2 x belong to the class L 2 (0, T; L 2 (fi)). 

By choosing T sufficiently small, it is possible to impose 



1 

2k 



dt\\h + h'VlP + 4 / dt (\\h' + Ac x c. 



T 



-'XX || 2 



1*4 + ^11!) 



< 1, 



which in turn forces sup T6 r 0T i H^/ijclli = 0^ an d hence the solution is unique. 
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